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I I < Abstract. - Transverse momentum spectra of charmed and beauty mesons produced in proton- 

i-^ ' proton and proton-antiproton collisions at high energies are analyzed within the modified quark- 
gluon string model (QGSM) including the internal motion of quarks in colliding hadrons. It is 

^ • shown that this approach can describe rather satisfactorily the experimental data at not large 

^ ' values of the transverse momentum where the NLO QCD calculation has a big uncertainty. We 

,.1^ , also show that using both the QGSM and the NLO QCD one can describe these data in a wide 

' region of transverse momenta and give some predictions for the future LHC experiments. 
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■ Introduction. — Various approaches of perturba- 
^^tive QCD including the next-to-leading order calculations 
('NLO QCD) have been apphed to construct distributions 
00 Qf quarks in a proton. The theoretical analysis of the lep- 
^^tx3n deep inelastic scattering (DIS) off protons and nuclei 
^provides rather realistic information on the distribution 
, " of Hght quarks like u,d,s in a proton. However, to find 
■ rN a reliable distribution of heavy quarks like c(c) and espc- 
\ oially 5(6) in a proton describing the experimental data on 
■ t'hc DIS is a non-trivial task. It is mainly due to small 
values of D and B meson yields in the DIS at existing 
energies. Even at the Tevatron energies the B- meson 
yield is not so large. At LHC energies the multiplicity of 
these mesons produced in pp collisions will be significantly 
larger. Therefore one can try to extract a new information 
on the distribution of these heavy quarks in a proton. In 
this paper we suggest to study the distribution of heavy 
quarks like c(c) and h{h) in a proton from the analysis of 
the future LHC experimental data. 

The multiple hadron production in hadron-nucleon col- 
lisions at high energies and large transfers is usually ana- 
lyzed within the hard parton scattering model (HPSM) 
suggested in [1, 2]. This model was apphed to the 
charmed meson production both in proton-proton and 
meson-proton interactions at high energies, see for exam- 
ple [3]. The HPSM is significantly improved by apply- 
ing the QCD parton approach implemented in the mod- 
ified minimal-subtraction renormalization and factoriza- 



tion scheme. The first calculation scheme is the so-called 
massive scheme or fixed-flavor-number scheme (FFNS) de- 
veloped in [4]- [7]. In this approach the number of ac- 
tive flavors in the initial state is limited to n/ = 4, e.g., 
u{u) , d{d) , s{s) and c(c) quarks being the initial partons, 
whereas the b{b) quark appears only in the final state. 
In this case the beauty quark is always treated as a heavy 
particle, not as a parton. In this scheme the mass of heavy 
quarks acts as a cutoff parameter for the initial-state and 
final-state collinear singularities and sets the scale for per- 
turbative calculations. Actually, the FFNS with Uf — A 
is limited to a rather small range of transverse momenta 
Pt of produced D oi B mesons less than the masses of c 
or b quarks. In this scheme the terms ^/pj are fully 
included. 

Another approach is the so-called zero-mass variable- 
flavor-number scheme (ZM-VFNS), see [8]- [10] and ref- 
erences therein. It is the conventional parton model ap- 
proach, the zero-mass parton approximation is also ap- 
plied to the b quark, although its mass is certainly much 
larger than the asymptotic scale parameter Aqcd- In this 
approach the b{b) quark is treated as an incoming parton 
originating from colliding hadrons. This approach can be 
used in the region of large transverse momenta of produced 
charmed or beauty mesons, e.g., at pt > rric^b- Within this 
scheme the terms of order ^/p^ can be neglected. Re- 
cently the experimental inclusive pt spectra of B mesons in 
pp collisions obtained by the CDFII Collaboration [12, 13] 
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at the Tevatron energy y/s = 1.96 TeV in the rapidity 
region — 1 < y < 1 have been described rather satisfac- 
torily within this ZM-VFNS approach in [11] at pt > 10 
GeV/c using the non-pcrturbative structure functions. In 
another kinematic region, e.g.. at 2.5 GeV/c < pt < 10 
GeV/c, the FFNS model allowed the CDFII data to be de- 
scribed without using fragmentation functions of b quarks 
to B mesons. Both these schemes have some uncertainties 
related to the renormalization parameters. 

In this paper we study the charmed and beauty me- 
son production within the QGSM [14] or the dual par- 
ton model (DPM) [15]. based on the 1/A^ expansion in 
QCD [16,17]. We show that this approach can be applied 
rather successfully at not very large values of pt . 

General formalism. — Let us analyze the _D-meson 
production in the pp and pp collisions within the QGSM 
including the transverse motion of quarks and diquarks 
in colliding protons [19]. As is known, the cylinder type 
graphs for the pp collision presented in fig. 1 make the 
main contribution to this process [14]. The left diagram 
of fig. 1, the so-called one-cyhnder graph, corresponds to 
the case where two colorless strings are formed between the 
quark/diquark (q/qq) and the diquark/quark (qq/q) in col- 
liding protons; then, after their breakup, qq pairs are cre- 
ated and fragmentated to a hadron, for example, D meson. 
The right diagram of fig. 1, the so-called multicylindcr 
graph, corresponds to creation of the same two colorless 
strings and many strings between sea quarks/ antiquarks 
q/q and sea antiquarks/quarks q/q in the colliding pro- 
tons. 
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Fig. 1: The one-cylinder graph (left diagram) and the multi- 
cylinder graph (right diagram) for the inclusive pp hX pro- 
cess. 
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Here r means the flavor of the valence (or sea) quark 
or diquark, fr"'\x' ^kt) is the quark distribution func- 
tion depending on the longitudinal momentum fraction 
x' and the transverse momentum kt in the n-Pomeron 
chain; Gr-^h{z, ktlPt) = zDr-,h{z, kt;pt), br^h{z,kt;pt) 
is the fragmentation function (FF) of a quark (antiquark) 
or diquark of flavor r into a hadron h {D meson in our 
case). We present the quark distributions and the FF in 
the factorized forms fr{x,kt) = fr{x)gr{kt), and, accord- 
ing to [22], Gr->h{z,kt]pt) = Gr^h{z)gT~*h(kt), where 
kt = Pt — ^kf. We take the quark distributions frix) 
and the FF Gr^h{z) obtained within the QGSM from 
[18,23,24], whereas their kt distributions are chosen in 
the form suggested in [21,22] 



The general form for the invariant inclusive hadron spec- 
trum within the QGSM is [18, 19] 



gAh) = (S„727r) exp(-B,( Jfc2 + ^1 - mo)) , (6) 



^d5^ = Vid^ ^ T.^-nis)Mx,Pt) , (1) 

where i?, p are the energy and the three-momentum of 
the produced hadron h in the laboratory system (l.s.) of 
colliding protons, E* , s are the energy of h and the square 
of the initial energy in the c.m.s of pp, x,pt are the Feyn- 
man variable and the transverse momentum of /i; (t„ is 
the cross section for production of the n-Pomeron chain 
(or 2n quark-antiquark strings) decaying into hadrons, cal- 
culated within the "eikonal approximation" [20] , the func- 
tion (j}n{x,pt) has the following form [19]: 



~gr^h{kt) = {BI/2T,) cM-Bc{^k1 + ml - mo)). (7) 

After the integration of eq. (5) over d^fct we have, ac- 
cording to [22], 

F^^\x±,pt;x^a) = h''\xi^2)Gr^h{z)In{z,pt) , (8) 

where z = x±/xi^2, In{z,Pt) = BI/{2'k{1 + 
B.mo)) exp{-B^{mDt - mo)), m|,j = Pt + = 
Be/ {I + npz^), p = Bc/Bq. The function B^ also can 
be presented in the equivalent form B^ ~ Bql[p + nz^). 
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where p = Bq/Bc- The differential cross section da /dpi 
for D mesons produced in pp collisions is written in the 
form [18] 



(9) 



The production of heavy mesons like D- and _B-mesons 
in proton-antiproton collisions at high energies is usually 
analyzed within the different schemes of QCD. To study 
these processes within the QGSM wc have to include at 
least one additional graph corresponding to the creation 
of three chains between quarks in the initial proton and 
antiquarks in the colliding antiproton, as is illustrated in 
fig. 2 (bottom diagram). 
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Tevatron energy, whereas experimental data on (7°^" are 
available only for the antiproton initial energy about 10 
GeV, see [26] and references therein. However, some theo- 
retical predictions, for example [27,28], show that asymp- 
totically (Jpp"' goes to about 2 — 4 mb. It corresponds 
to (jj ~ '^pp^ l<^]°p < 0.1 at the Tevatron energy. Note 
that in addition to the graph of fig 2c there can be di- 
agrams consisting of these three chains and multicilyn- 
der chains between sea quarks and antiquarks. However, 
as our estimations show, their contribution to the in- 
clusive spectrum is much smaller than the contribution 
from the three-chain graph (fig 2c). Therefore we ne- 
glect it. The form for the function (t)^{x,pt) is similar 
to 4>n{x^pt) entering into (3) by replacing Fq^\x-,pt', X2), 



Fq'^\OiPt) to F^q'{x-,pt;x2), F^q'{Q,pt) respectively, and 
replacing F^q\x-,pt;x2), F!^q\o,pt) to F^''\x-,pt; X2) 
and Fq"'\o,pt) respectively. The additional term (f>{x,pt) 
in (11) has the following form 



(")/ 



4>ix,Pt) = 3Fq^{x+,Pt)Fq^{X-,Pt)/Fq{0,Pt) 



(11) 
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Fig. 2: The one-cylinder graph (left diagram) and the multi- 
cylinder graph (right diagram), and the three-chains graph 
(down diagram) for the pp hX inclusive process. 

The left and right diagrams in fig. 2 are similar to the 
one-cylinder and multicylinder diagrams for the pp colli- 
sion in fig. 1 with a following difference. In the pp col- 
lision two colorless strings between quark/diquark {q/ qq) 
in the initial proton and antiquark/antidiquark (q/qq) are 
created. Many quark-antiquark {q — q) strings for pp col- 
lision (fig. 2, right diagram) are the same as for the pp 
collision (fig. 1, right diagram). Therefore, the invariant 
inclusive spectrum of hadrons produced in the pp collision 
calculated within the QGSM has the following form: 



= ai(s)((l - co)0f (x,ft) + iucj,ix,pt)) + (10) 



d3p 



^an{s)(l)PP{x,pt) 



n=1 

where 1 — w is the probability of contribution of the 
cut one-cylinder (one-Pomeron exchange) and cut multi- 
cylinder (multiPomeron exchanges) graphs (the left and 
right diagrams in fig. 2), whereas lo is the probability of 
the contribution of the three-chain diagram to the inclu- 
sive spectrum. The value of lo can be estimated as the 
ratio of the pp annihilation cross section a^l^" to the to- 
tal pp cross section cr*^*. The cross section cr*^* is well 
known in the wide range of the initial energies to the 
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Results and discussion. — To illustrate our ap- 
proach we present in figs. 3,4 the inclusive spectrum 
da /dpi of D°-mesons produced in the reaction D^X 
at -^s ~ 27.4 GcV as a function of pi given by eq.(9). 
One can see from fig. 3 that the use of different values 
for the intercept a*(0) leads mainly to the shift of the 
theoretical lines along the y axis. figs. 3,4 show a satis- 
factory description of the experimental data obtained by 
the NA27 Collaboration [25] at both values for the linear 
VP-trajectory when a\s,{0) ~ —2.18 and for the nonlinear 
one when a^(0) = 0. Unfortunately, the experimental 
data are very poor because of big error bars; therefore, we 
cannot get new information on the Vf-trajectory. 

The inclusive pt spectra of 13° and B~^ mesons produced 
in the pp collision at the Tevatron energy y/s = 1.96 
TcV are presented in figs. 5,6. The hatched regions in 
figs. 5,6 show the calculations within the NLO QCD in- 
cluding uncertainties [30]. The dashed hues (w = 0) and 
dash-dotted curves (w = 0.1) in figs. 5,6 correspond to our 
best fit obtained within the QGSM-I calculations using 
the parameter values Be = 0.65(GeV/c)^^ for D mesons 
(fig. 5) and B^ = 0.55(GeV/c)~^ for B mesons (fig. 6), 
and p = Bc/Bq = 3.1 [22], e.g., Bq ~ 0.2(GeV/c)-i. 
However, the value used for the slope Bq of the quark dis- 
tribution as a function of kt is too small. Therefore, we 
also calculated these pt-spcctra taking the more realistic 
values Bq = 4.5(GeV/c)~^ at the Tevatron energy and 
Bq = 4(GeV/c)^^ at the LHC energy that correspond ap- 
proximately to < kt >— 0.45 GeV/c and < kt >— 0.5 
GeV/c respectively. This calculation (QGSM-II) is shown 
by the solid lines in figs. 5,6 and in figs. 7,8. Accord- 
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Fig. 3: The inclusive spectrum for mesons produced in the 
pp collision at ^/s — 27 A GeV as a function of pi , see details 
in [22]. 



Fig. 4: The inclusive spectrum for Z)° mesons produced in the 
pp collision at ^/s = 27.4 GeV as a function of x, see details 
in [22]. 



ing to the experimental data, the mean transverse mo- 
mentum of hadrons produced in hadron collisions slowly 
increases as the energy increases. Therefore, the inter- 
nal transverse momentum of the quark in the proton can 
also slow increase. In fact, we have only one free param- 
eter p= 1/ p = Bq/Bc which is experimentally unknown, 
whereas Bq is directly related to the mean transverse mo- 
mentum of the quark in the proton or antiproton which 
is more or less known experimentally. To describe the ex- 
perimental data on the pt-spectra in the pt region where 
the NLO QCD calculation has a big uncertainty we chose 
p = 7 both for the Tevatron and the LHC energies. When 
Be < I (GeV/c)~^, eq.(7) can be approximately pre- 
sented in the form gr^hikt) — a^/(a^ -I- kf) at < 2mD, 
where a = ^ImDjBf^. This form is similar to the form for 
the FF of heavy quarks obtained within the perturbative 
QCD, see for example [29]. Note that the function In{z,pt) 
in eq.(8) was obtained in [19,22] on the assumption of 
the consequent sharing of the transverse momentum pt in 
the proton (antiproton) between n-Pomeron chains. It al- 
lowed us to describe rather satisfactorily the experimental 
data on the inclusive pt spectra of charmed and beauty 
mesons produced in pp collisions at moderate values of 
the transverse momentum pt < 10 GeV/c. It is illus- 
trated in figs. 5,6. We also use A ~ 2a'^.^^.^(0) < pf >, 
"ri-(B>-)(0) - 0.5 (GeV/c)-2 is the slope of the D* or B* 



Regge trajectory, < p^ >~ 5 (GeV/c)^ is the mean trans- 
verse momentum squared of the D meson or B meson that 
was found from the CDFII experimental data. Note that 
the our calculation showed that the contribution of the 
three-chain graph (fig. 2c) is very small at the Tevatron 
energy. It is due to small values of the pp annihilation 
cross section at very high energies [27,28]. 

The predictions for inclusive pt spectra of and 5+ 
mesons produced in the pp collision at LHC energies 
and the NLO QCD calculation for the produced charmed 
quarks [31] are presented in figs. 7,8. 

The solid lines correspond to our calculations within the 
QGSM-II, whereas the hatched regions show the calcula- 
tions within the NLO QCD including uncertainties [31]. A 
big difference between the QGSM and NLO QCD calcula- 
tions at Pt > 10 GeV/c for D and B mesons can be due to 
the following. First, the NLO QCD calculation [31] does 
not include the hadronization of quarks to heavy mesons, 
whereas the QGSM calculation includes it. Second, we 
do not include the contribution of gluons and their hard 
scatterings off quarks and gluons which can be sizable at 
large values of pt . 

Conclusion. — We have shown that the modified 
QGSM including the intrinsic longitudinal and transverse 
motion of quarks (antiquarks) and diquarks in colliding 
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Fig. 5: The inclusive pt-spectrum for D mesons produced 
in the pp collision at the Tevatron energy ^/s — 1.96 TeV 
obtained within the QGSM (the solid and dashed lines) and 
within the NLO QCD [30] (the hatched regions); QGSMT: 
Bg ~ 0.18(Gev/c)-\ p ~ 0.32; QGSMTI: = 4.5(GeV/c)~\ 



protons allowed us to describe rather satisfactorily the ex- 
isting experimental data on inclusive spectra of D mesons 
produced in pp collisions and to make some predictions for 
similar spectra at LHC energies. To verify whether these 
predictions can be reliable or not we apply the QGSM to 
the analysis of charmed and beauty meson production in 
proton- antiproton collisions at Tevatron energies including 
graphs like those in fig. 2c corresponding to annihilation of 
quarks and antiquarks in colliding p and p, and production 
of D-mesons. 

We got a satisfactory QGSM-II description {pt < 10 
GeV/c) of the experimental data on pt spectra of and 
mesons produced in the pp collisions which were ob- 
tained by the CDFII Collaboration at the Tevatron [30]. 
At larger values of pt the calculations within the NLO of 
QCD [4, 8] result in a better description of these data. 
It can be due to the contribution of gluons inside the 
colliding proton and antiproton which can interact with 
other gluons and quarks (antiquarks) and fragmentate to 
charmed mesons. This effect is not taken into account in 
the presented QGSM. Therefore, the QGSM including the 
internal transverse momenta of partons in the proton (let 



Fig. 6: The inclusive pt-spectrum for B mesons produced 
in the pp collision at the Tevatron energy = 1.96 TeV 
obtained within the QGSM (the solid and dashed lines) and 
within the NLO QCD [30] (the hatched regions); QGSM-L 
Bq ~ 0.21(GeV/c)-\ p ~ 0.32; QGSM-II is the same as in 
Fig (5). 

us call it the "soft QCD") allows us to describe the in- 
clusive Pt spectra of heavy mesons produced in pp and pp 
collisions at high energies at not very large values of pt- 
To describe these spectra and make some predictions for 
the future LHC experiments in a wide region of transverse 
momenta one can combine the "soft QCD" at small values 
of with the NLO QCD at large pt- 

We found that the pt spectra of D and B mesons calcu- 
lated within the QGSM are almost insensitive to the form 
of the sea c(c) and b{b) quark distributions in colliding 
protons/antiprotons. To find a new information on it we 
intend to study the charm and beauty hadron production 
in pp collisions at small scattering angles. 
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Fig. 7: The inclusive spectrum for mesons produced in 
the pp colhsion at the LHC energy -^s — 14 TeV obtained 
within the QGSM for charmed mesons and the NLO QCD for 
c quarks [31]; QGSM-I corresponds to the same QGSM-I as in 
Fig 5; QGSM-II; = 4. (GeV/c)-\ p = 7. 
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